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Another book recommendation 


Latitude 

How American Astronomers Solved the Mystery of Variation 


Bill Carter and Merri Sue Carter 


By Bill Carter and Merri Sue Carter, 
Naval Institute Press, Annapolis, 
Maryland 2002 



















Review 


Latitude and Longitude can uniquely and 
meaningfully describe where we are on 
the earth. 

We can also express positions on a 
sphere using 3-D Cartesian coordinates 
[ X;Y;Z] using simple geometric 
relationships. 





























Textbook error 


See page 23. 

Author misplaces the decimal point where he 
converts DMS to Decimal 

- 31.315278 should be 3.1315278 
The answer to problem 2.2 is correct. 



Homework Answer 


• Problem: Use the position of station BLUCHER 
to determine the distance on a spherical earth 
(with radius 6,378,000 m) from the equator. 

- BLUCHER: 27-42-52.08857N 



Both methods yield: 3,085,094 m 























What about using INVERSE? 


Output from INVERSE 


Ellipsoid : GRS80 / WGS84 
Equatorial axis, a 

Polar axis, ta 

Inverse flattening, 1/f = 


(NAD83) 

6378137.0000 
6356752.3141 
298.25722210088 


First Station : equator 


LAT =00 0.00000 North 

LON = 97 19 44.31265 West 


Second Station : BLUCHER 


LAT = 27 42 52.08857 North 

LON = 97 19 44.31265 Uest 


Forward azimuth FAZ 
Back azimuth BAZ 
Ellipsoidal distance S 


0 0 0.0000 From North 
180 0 0.0000 From North 
3066800.0198 m 



3,066,800 meters 
10,061,660 feet (US Survey) 
1,905.6 miles (statute) 


Our result using a radius of 
6,378,000 meters is 3,085,094 m 

A difference of 18,294 m! 

Why?? 










Feet are Feet? 


Conversions from meters to feet (and 
inverse) are complicated by two units of 
feet. 

- U.S. Survey foot = 0.30480061... meters 

• 1200/3937 meters (exactly) 

- International foot = 0.3048 meters (exactly) 

• 2.54 cm = 1 inch 



DMS <-> Radian 

To convert degrees to radians 

- Convert DD MM SS.sssss to decimal 

• Deg + m in/60 + sec/3600 

- Convert decimal degrees to radians 

• Multiply by pi/180 

To convert radians to decimal 

-decDeg = Radian value * 180/pi 

- Deg = floor(decDeg) 

- Min = floor((decDeg-Deg)*60) 

- Sec = decDeg*3600-(Deg*60)-(Min*3600) 



The meter 


There were great difficulties in commerce due to 
varying length (and other) units. 

The French Academy of Science was charged 
with standardizing the measurement unit. 

Original proposal was to use the period of a 
pendulum. 

Instead, in 1790 the Academy recommended 
that a meter unit be based on one-millionth of 
the distance from the Equator to the North Pole. 



How well did they do? 


Output from INVERSE 


Output from INVERSE 


Ellipsoid : Clarke 

1866 

(NAD27) 

Equatorial axis. 

a = 

6378206.4000 

Polar axis. 

b 

6356583.8000 

Inverse flattening. 

1/f = 

294.97869821380 


Ellipsoid : GRS80 / 

WGS84 

(NAD83) 

Equatorial axis. 

a = 

6378137.0000 

Polar axis. 

b 

6356752.3141 

Inverse flattening. 

1/f = 

298.25722210088 


First Station : 


First Station : 


LAT =00 0.00000 North 

LON = 100 0 0.00000 West 

Second Station : 


LAT = 90 0 0.00000 North 

LON = 100 0 0.00000 West 


Forward azimuth 

FAZ = 0 

0 

0.0000 

From 

North 

Back azimuth 

BAZ = 180 

0 

0.0000 

From 

North 

Ellipsoidal distance 

S = 10001888.0430 

HI 



LAT =00 0.00000 North 

LON = 100 0 0.00000 West 

Second Station : 


LAT = 90 0 0.00000 North 
LON = 100 0 0.00000 West 


Forward azimuth 

FAZ = 0 

0 

0.0000 

From 

North 

Back azimuth 

BAZ = 180 

0 

0.0000 

From 

North 

Ellipsoidal distance 

S = 10001965.7292 

HI 


















Evolution of the meter 


The original measurement was in error due to 
unknown magnitude of Earth’s flattening. 

- The unit was transferred to a platinum-iridium alloy 
bar kept in Paris (1874) 

- The unit was updated in (1889) to a bar composed of 
90% platinum 

In 1960 a new definition was adopted that was 
based on krypton-86 radiation wavelength. 

Meter is the length of the path traveled by light in 
a vacuum during the time interval of: 

- 299 792 458 s-1 ( 299 792 458 meters per sec) 



Distance on a sphere 



Compute radius of parallel circle by solving right triangle. 





















Subtract longitudes to get angle, s = Rp * angle (in radians) 
For this problem s = 179,237 meters 





Spherical Triangles 


Used in great circle navigation. 

Sides and angle are measured using arc 
measures 

Located on the surface of the sphere with 
sides formed by great circle arcs. 

- N.B. great circles are planes through the 
center of the Earth 

The shortest distance between points. 

- Not exactly 



Spherical trigonometry 









More damned definitions 


Normal section is a plane that contains the 
normal to the sphere at the occupied point and 
another point of interest. 

Horizontal angle is measured between two 
normal sections with respect to the instrument 
location. 

Azimuth is measured from the normal section 
containing the N pole clockwise to the normal 
section containing the other point 

All normal sections on the sphere intersect 
the sphere along great circle arcs. 



Spherical Triangle 


A X = Jricluded Angle of 
Spherical Triangle 


Rotation Ajus 


Flans 
of 
Psmllel 


Trace 


Great 

Circle 



Co-ktitodle = Sides of 
Spherical Triangle 










Problem in Text 







Distance Calculations 

Determine the co-latitudes (90d - latitude) at 
points A and B. 

- These are the lengths C to B (side a) and C to A (side 
b). 

Compute the difference in Longitudes 

- This is the angle at C 

As Law of Sines is ambiguous for angles in 
excess of 90d, we use Law of Cosines to solve 
for distance side c 

- Cos(C) = cos(A)cos(B)+sin(A)sin(B)cos(C). 

Distance = r * C (N.B. radius of great circle same 
as the sphere itself. 



We can also calculate Azimuth 


Use cotangent formulas and the results 
from the spherical triangle computation. 
-tanA = sinC/( (sin(b)/tan(a))-(cos(b)cos(C)) 

Note that we must correctly account for 
the quadrant. 

Note as well that forward and reverse 
azimuths are not exactly 180 d different. 



Other spherical Earth 
characteristics 

All meridians converge at poles. 

Azimuths of lines measured from one end to not 
equal values measured from the other end. 

- Effect is especially pronounced on long E-W lines. 

Can be approximated as a function of the E-W 
distance, mean latitude and spherical radius. 


r p dXany/ 
R 





Spherical excess 


The summation of all spherical angles 
exceed 180 degrees. 

It is proportional to the area of the 
spherical triangle. 


be sin A 
2R 2 sin 1" 




|aP20 0 6 
1356 
29 


## 

+ 

+ 

+ 

+ 

+ 

+ + 
+ + 
+ + 
+ + 
+ + 
%c 
%c 
%f 
%f 
%i 
%i 
/* 


1100 0.00000000 96 ORBIT IGbOO HLM IGS 

0.00000000 900.00000000 53736 0.0000000000000 

1 2 3 4 5 6 7 8 9 10 11 13 14 15 16 17 18 
19 20 21 22 23 24 25 26 27 28 29 30 0 0 0 0 0 

0 


0 

0 

0 

3 

4 
0 
0 
0 


0 

0 

0 

3 

3 

0 

0 

0 


0 

0 

0 

4 

3 

0 

0 

0 


0 

0 

0 

3 

3 

0 

0 

0 


0 

0 

0 

3 

3 

0 

0 

0 


0 

0 

0 

3 

3 

0 

0 

0 


0 

0 

0 

3 

4 
0 
0 
0 


0 

0 

0 

3 

3 

0 

0 

0 


0 

0 

0 

3 

3 

0 

0 

0 


0 

0 

0 

3 

4 
0 
0 
0 


0 

0 

0 

3 

3 

0 

0 

0 


0 

0 

0 

3 

3 

0 

0 

0 


0 

0 

3 

0 

0 

0 

0 


0 

0 

0 

3 

0 

0 

0 

0 


0 

0 

0 

4 

0 

0 

0 

0 


0 

0 

0 

3 

0 

0 

0 

0 


0 

0 

0 

3 

0 

0 

0 

0 


cc cc 
cc cc 
0 . 00 
0 . 00 
0 
0 

FINAL 


ccc ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc 
ccc ccc cccc cccc cccc cccc ccccc ccccc ccccc ccccc 
00000 0.000000000 0.00000000000 0.000000000000000 
00000 0.000000000 0.00000000000 0.000000000000000 
0000000 0 
0000000 0 
ORBIT COMBINATION FROM WEIGHTED AVERAGE OF: 


/* cod emr esa gfz jpl mit ngs sio 

/* REFERENCED TO IGS TIME (IGST) AND TO WEIGHTED 


/* CLK ANT Z-OFFSET 
* 2006 1 1 0 0 
f P 1 -15156. 118389 

P 2 20133.793284 

P 3 -11366.447142 
P 4 26109.950844 

P 5 -11013.232261 
P 6 -10108.117501 
P 7 16494.168125 

P 8 26002.033551 

P 9 -872.962985 

P 10 5783.888427 

P 11 2597.964722 


(M) : II/IIA 1, 
0.00000000 


023; HR 0.000 


18164.U21963 -11829.8875881 


-12517.384408 11531.263733 

22473.599105 7788.691364 

-5714.160584 1216.379808 

-21731.530250 -10967.998519 
-15606.106217 19125.583219 

-8125.986115 -18697.734751 
50.071986 6220.016726 

-15641.294915 -21956.384902 
-15080.250565 20946.780855 

16510.768220 -20675.855544 


MEAN POLE: 


28.085474 
-23.264373 
63.023929 
103.694692 
447.812380 
170.870195 
459.605468 
-51.940408 
1.771629 
76.347737 
295.573462 






£. 10 
teqc 2000Feb29 
MSWinNT 4.0|Pent Pro|bcc32 


OBSERVATION DATA 


BIT 2 OF LLI 
-Unknown- 
-Unknown- 
-Unknown- 
-Unknown- 

-726296.8700 
0.0000 
1 1 
7 LI 


G (GPS) 

20050912 16:17 
0 | MSWin95/98/NT|4 8 6/DX+ 
FLAGS DATA COLLECTED UNDER A/S CONDITION 


-Unknown- 
ASHTECH UZ-12 
-Unknown- 

-5598342.1900 2958518 

0.0000 0 


L2 Cl P2 PI 
SNR is mapped to RINEX snr flag value 

LI: 1 -> 1; 90 -> 5; 210 -> 9 

L2: 1 -> 1; 150 -> 5; 250 -> 9 


ZE21 

6100 

0000 

Dl 

[1-9] 


D2 


2005 


14 


13 


37 


15.0000000 


GPS 


RINEX VERSION / TYPE 
4OUTCPGM / RUN BY / DATE 
COMMENT 
COMMENT 
MARKER NAME 
OBSERVER / AGENCY 
REC # / TYPE / VERS 
ANT # / TYPE 
APPROX POSITION XYZ 
ANTENNA: DELTA H/E/N 
WAVELENGTH FACT Ll/2 
# / TYPES OF OBSERV 
COMMENT 
COMMENT 
COMMENT 

TIME OF FIRST OBS 
END OF HEADER 


05 9 14 13 37 15. 

0000000 0 6G 

5G30G10G18G26G 

2 


138416.07353 

-2389.228 

79208.19652 

-1861.736 

21658215.946 

21658212.382 

21658216.095 

128476.45353 

-2125.583 

92042.54252 

-1656.298 

20081603.431 

20081600.703 

20081604.871 

179425.29253 

-2613.179 

29628.04052 

-2036.243 

19575418.616 

19575412.115 

19575414.890 

-154209.07053 

2887.546 

-111215.70852 

2250.036 

20094561.440 

20094559.983 

20094564.353 

-137048.39353 

2422.747 

-93602.12652 

1887.855 

17903365.528 

17903361.718 

17903367.264 

113885.28653 

-1716.686 

82360.75052 

-1337.677 

20632766.604 

20632762.002 

20632768.438 



